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1.0  SUMMARY 


The  overarching  goal  of  this  project  was  initially  to  study  how  obesity  modifies  the 
effects  of  genetic  variants  within  previously  identified  genes  for  asthma  using  previously 
collected  genotype  and  phenotype  data  from  a  number  of  existing  genetic  epidemiological 
studies.  Through  the  course  of  this  project,  we  made  some  interesting  observations  that  led  us  to 
expand  the  scope  of  our  search  to  take  a  genome-wide  approach  rather  than  a  candidate  gene 
approach.  Further,  after  identifying  several  genome-wide  significant  interactions  between  body 
mass  index  (BMI)  and  single  nucleotide  polymorphisms  (SNPs)  associated  with  asthma,  we 
looked  at  these  SNPs  in  more  detail  to  more  fully  understand  the  relationship  between  the  SNPs 
and  obesity  and  asthma.  We  made  an  interesting  observation  that  the  SNPs  appeared  to  have  an 
effect  on  BMI  among  asthmatics  but  not  non-asthmatics.  This  observation  was  made  among 
U.S.  whites  in  the  Multi-Ethnic  Study  of  Atherosclerosis  (MESA)  dataset  but  not  for  other  ethnic 
groups  in  MESA  and  Women’s  Health  Initiative  (WHI).  None  of  the  top  SNPs  overlapped 
among  any  of  the  ethnic  groups.  The  results  of  this  analysis  altered  the  direction  of  this  project 
so  that  we  can  better  understand  the  relationship  between  these  SNPs,  BMI,  and  asthma. 

We  started  by  looking  at  BMI  as  an  outcome  among  the  U.S.  whites  from  MESA.  To 
check  if  interactions  between  significant  SNPs  and  asthma  status  in  MESA  are  replicable,  the 
Framingham  Heart  Study  (FHS)  was  used  as  an  independent  dataset  to  replicate  the  top  SNPs  in 
MESA  due  to  the  availability  of  subjects  of  similar  age  and  ethnicity.  There  are  three  generations 
in  the  Framingham  study.  The  second  generation,  the  offspring  sub-study,  was  chosen  for 
analysis  because  its  participants  have  similar  age  (60-70  years  old)  as  MESA. 

We  performed  a  genome-wide  gene  by  environment  (asthma)  interaction  analysis  for  the 
outcome  of  BMI  in  MESA  (N=2474  Caucasians,  257  asthmatics)  and  replicated  findings  in  the 
FHS  offspring  cohort  (N=1408  Caucasians,  382  asthmatics).  The  replicable  tagging  SNP 
rs2107212  was  further  examined  in  stratified  analyses.  Seven  SNPs  clustered  in  1 7q2 1 .2  were 
identified  to  be  associated  with  higher  BMI  among  asthmatics  (interaction  p  <  5  x  10"7  in  MESA 
and  p  <  0.05  in  FHS).  In  both  MESA  and  FHS  asthmatics,  subjects  carrying  the  A  allele  on 
rs2107212  had  significantly  higher  odds  of  obesity  than  non-carriers,  which  was  not  the  case  for 
non-asthmatics.  We  further  examined  BMI  change  subsequent  to  asthma  diagnosis  over  a  period 
of  26  years  in  FHS  and  demonstrated  greater  BMI  increase  among  asthmatics  compared  to  non¬ 
asthmatics.  Asthmatics  carrying  the  A  allele  at  rs2107212  had  significantly  greater  net  BMI 
increase  over  the  26-year  period  compared  to  non-asthmatics.  In  this  study,  we  found  that 
common  genetic  variants  on  1 7q2 1 .2  are  associated  with  post-asthma  BMI  increase  among 
Caucasians. 

This  work  is  detailed  in  a  published  peer-reviewed  paper:  Wang  L,  Murk  W,  DeWan  AT. 
Genome-wide  gene  by  environment  interaction  analysis  identifies  common  SNPs  at  1 7q2 1 .2  that 
are  associated  with  increased  body  mass  index  only  among  asthmatics.  PLoS  One.  2015; 
10(12):e01441 14.  doi:  10.1371/journal.pone.01441 14.  eCollection  2015. 

We  further  explored  the  relationship  between  a  different  type  of  genetic  variant,  copy 
number  variants  (CNVs),  which  exist  as  deletions  and  duplications,  asthma,  and  obesity.  Using 
the  WHI  dataset,  we  analyzed  the  African  American  subjects  (due  to  power  constraints)  for 
interactions  between  CNVs  and  BMI  that  were  associated  with  asthma.  The  most  significant 
deletion  interacting  with  BMI  associated  with  asthma  is  on  chromosome  7  (p-value  =  0.0032). 
The  most  significant  duplication  interacting  with  BMI  in  the  outcome  of  asthma  is  on 
chromosome  10  (p-value  =  0.0019). 
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2.0  INTRODUCTION 


Asthma  and  obesity  are  two  rapidly  growing  public  health  issues,  and  the  comorbidity  of 
these  conditions  poses  an  enormous  burden  on  asthma  control  as  well  as  quality  of  life  [1].  The 
simultaneously  increasing  prevalence  of  both  asthma  and  obesity  suggests  a  potential  intrinsic 
link  between  these  two  chronic  disorders  [2,3],  and  a  recent  longitudinal  study  identified  asthma 
as  a  risk  factor  for  subsequent  obesity  [4],  However,  the  underlying  factors  that  contribute  to  this 
relationship  have  remained  largely  unknown. 

Genetic  factors  play  an  essential  role  in  both  asthma  and  obesity  and  are  believed  to  be 
predominantly  responsible  for  the  comorbidity  of  the  two  conditions  [5-7].  A  genomic  inversion 
in  16pl  1.2  was  identified  to  be  protective  against  the  joint  occurrence  of  asthma  and  obesity  in 
adults  of  European  descent  [8],  This  provided  genetic  evidence  of  asthma-obesity  co-occurrence. 
In  terms  of  single  nucleotide  polymorphisms  (SNPs),  a  candidate  gene  association  study  was 
conducted  to  identify  shared  genetic  variants  between  childhood  asthma  and  obesity,  but  no  SNP 
was  associated  with  both  phenotypes  among  previously  identified  asthma  and  body  mass  index 
(BMI)  genes  [9].  This  implied  that  SNPs  and  copy  number  variants  (CNVs)  underlying  the 
comorbidity  of  asthma  and  obesity  may  exist  in  novel  loci.  A  genome-wide  analysis  study 
(GWAS)  was  also  conducted  to  identify  genetic  variants  associated  with  BMI  among  23,000 
asthmatics  [10].  An  SNP  in  DENND1B  was  identified  in  asthmatic  children  in  the  discovery 
dataset  but  was  not  replicable.  These  two  SNP  studies  primarily  targeted  childhood  asthma.  A 
potential  explanation  of  the  non-significant  findings  may  be  that  the  association  between  BMI 
and  asthma  is  not  as  strong  as  in  adults,  since  no  significant  BMI  difference  between  asthmatics 
and  non-asthmatics  was  observed  at  baseline  in  the  GWAS  study. 

3.0  METHODS 

3.1  Subjects  and  Phenotype  Definition 

The  Multi-Ethnic  Study  of  Atherosclerosis  (MESA)  and  its  ancillary  study,  MESA  Air, 
were  used  for  SNP  analysis  [11],  MESA  is  a  population-based  study  focusing  on  characteristics 
and  risk  factors  of  subclinical  cardiovascular  disease.  The  study  comprises  6,814  men  and 
women  aged  45-84  who  were  free  of  clinical  cardiovascular  diseases,  recruited  through  six  field 
centers  across  the  United  States.  The  screening  exam  (examl)  took  place  in  2000  and  was 
followed  by  four  examinations  (exams  2-5)  in  2002,  2004,  2005,  and  2010,  respectively.  Data 
from  exams  1-4  were  available  at  the  start  of  the  project  and  were  used  for  analysis.  Height  and 
weight  were  measured  at  every  visit,  and  the  BMI  value  [BMI  =  weight  (kg)/height  (m)2]  from 
the  screening  exam  was  used  in  the  present  analysis.  Asthmatics  were  defined  as  those  who 
reported  doctor-diagnosed  asthma  in  the  first  exam;  non-asthmatics  were  defined  as  those  who 
never  reported  doctor-diagnosed  asthma  or  asthma  medication  use  in  all  four  exams.  We 
included  all  participants  in  the  first  examination  who  identified  themselves  as  Caucasian 
(n=2527).  Subjects  with  missing  information  on  asthma  phenotype  (n=10),  inconsistent 
phenotype  infonnation  (i.e.,  those  who  used  asthma  medication  but  did  not  report  doctor- 
diagnosed  asthma,  n=134),  or  incident  asthma  after  first  examination  (i.e.,  those  who  first 
reported  asthma  in  exams  2-4,  n=34)  were  excluded.  MESA  Air  added  253  new  participants,  and 
after  exclusion  of  subjects  with  missing  infonnation,  we  included  249  additional  subjects.  Of  the 
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remaining  subjects  from  both  MESA  and  MESA  Air  (n=2598),  2588  were  genotyped  as  part  of 
the  MESA  SNP  Health  Association  Resource  (SHARe)  and  were  included  for  further  analysis. 

Data  from  the  Framingham  Heart  Study  (FHS)  were  used  for  replication  of  the  top  SNPs 
from  the  discovery  stage  [12].  The  objective  of  FHS  is  to  identify  common  factors  or 
characteristics  that  contribute  to  cardiovascular  disease.  It  was  established  in  1948  and  has 
developed  as  a  prospective,  community-based,  three-generation  study.  We  selected  participants 
in  the  offspring  cohort  for  analysis,  since  this  population  had  a  similar  age  distribution  as  the 
MESA  population.  There  were  eight  exams  from  1971  to  2005,  and  height  and  weight 
information  was  collected  at  each  visit.  BMI,  calculated  based  on  the  height  and  weight 
information  in  the  eighth  examination  (BMI  =  weight  (lb)/height  (in)2  x  703),  was  used  in  the 
present  analysis.  Except  for  the  first  exam,  self-reported  asthma  and  clinical  diagnostic 
impression  of  asthma  (GDI  asthma)  questions  were  asked  at  each  visit.  Asthmatics  were  defined 
as  those  who  ever  had  self-reported  asthma  or  wheezing  in  any  exams.  Non-asthmatics  were 
defined  as  those  who  never  had  self-reported  asthma  and  wheezing  or  GDI  asthma.  We  included 
subjects  with  non-missing  height  and  weight  information  in  the  eighth  examination  (n=2852). 
Subjects  with  GDI  asthma  but  not  self-reported  asthma  (n=13),  without  genotype  data  (n=178), 
or  not  confirmed  as  Caucasians  (n=3 1 1)  were  removed,  resulting  in  a  total  of  2350  subjects 
included  in  further  analysis. 

For  the  CNV  analysis,  we  used  genotype  and  phenotype  data  from  the  Women’s  Health 
Initiative  (WHI).  WHI  is  a  long-term  national  health  study  including  three  clinical  trials  and  one 
observational  study.  The  goal  of  the  observational  study  is  to  identify  risk  factors  associated  with 
specific  health  or  disease  outcomes.  Subjects  enrolled  in  the  WHI  were  postmenopausal  women 
between  50  and  79  years  old,  likely  to  live  in  the  area  for  3  years  and  have  no  condition 
predicting  survival  <3  years.  The  enrollment  was  completed  in  1998  and  participants  were 
followed  annually  for  8  to  12  years.  Genotyping  was  performed  using  Affymetrix  array  on 
African  Americans  (n=8218)  and  Hispanics  (n=3482).  CNVs  were  detected  through  genotype 
infonnation  for  African  Americans  due  to  power  constraints  for  the  Hispanic  subjects. 

3.2  Ethical  Statement 

Data  for  both  MESA,  FHS,  and  WHI  were  obtained  from  the  database  of  Genotypes  and 
Phenotypes  (MESA  accession  number:  phs000209.vl2.p3;  FHS  accession  number: 
phs000007.v23.p8;  WHI  accession  number:  phs000200.vl0.p3).  This  study  was  approved  by 
Yale  University  Human  Investigation  Committee.  Patient  records/infonnation  were  anonymized 
and  de-identified  prior  to  analysis. 

3.3  Genotyping  and  Quality  Control 

In  MESA,  genomic  DNA  was  isolated  from  peripheral  blood  samples  and  genotyped  on 
the  Affymetrix  Genome-Wide  Human  SNP  Array  6.0  chip  containing  909,622  SNPs.  Quality 
control  procedures  and  subsequent  genetic  analyses  were  performed  using  PLINK  1.07  [13]. 
Nineteen  subjects  were  removed  because  their  overall  call  rate  was  <98%,  and  six  subjects  were 
removed  because  the  reported  sex  did  not  match  the  genotyped  sex  as  called  in  PLINK.  SNPs 
were  excluded  if  they  met  any  one  of  the  following  criteria:  1)  not  mapped  to  an  autosomal 
chromosome  (n=37,380);  2)  call  rate  <98%  (n=2 1,669);  3)  the  genotype  distributions  of  the  SNP 
deviated  from  those  expected  by  Hardy-Weinberg  equilibrium  at  p  <  lxlO'7  (n=  1,729); 
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4)  monomorphic  or  minor  allele  frequency  (MAF)  <1%  (n=126,951).  A  subset  of  linkage 
disequilibrium  (LD)-pruned  SNPs  was  generated  using  pairwise  r2>0.5  in  a  window  of  50  SNPs 
and  shifted  by  5  SNPs  at  each  step  for  genetic  quality  control  procedures,  including  identity-by¬ 
descent  and  principal  component  analyses.  The  pairwise  identity-by-descent  matrix  was 
calculated  to  assess  cryptic  relatedness.  Thirty-three  pairs  of  subjects  were  identified  that  had  a 
pi-hat  value  >0.2.  One  subject  from  each  pair  was  randomly  removed.  Principal  component 
analysis  was  conducted  using  default  parameters  in  EIGENSTRAT  3.0  [14]  to  address 
population  stratification.  Fifty-six  outliers  were  identified  after  five  iterations  and  excluded  from 
further  analysis.  In  total,  2474  subjects  (n=257  asthmatics,  n=2217  non-asthmatics)  with  721,893 
SNPs  were  included  in  the  genome-wide  interaction  association  analysis  (discovery  stage). 

In  the  Framingham  SHARe,  a  sub-study  of  FHS,  subjects  were  genotyped  on  the 
Affymetrix  500K  mapping  array  plus  Affymetrix  50K  supplemental  array,  with  500,568  SNPs  in 
total.  The  same  quality  control  procedures  and  criteria  as  used  in  MESA  were  applied  to  the  FHS 
dataset.  Subjects  with  call  rate  <98%  (n=277)  and  subjects  with  unmatched  sex  check  (n=8) 
were  removed.  SNPs  that  were  not  mapped  to  an  autosomal  chromosome  (n=  12,422),  had  a  call 
rate  <98%  (n=79,984),  failed  Hardy- Weinberg  equilibrium  (n=l,776),  or  had  MAF  <1%  or  were 
monomorphic  (n=57,689)  were  excluded.  LD-pruned  SNPs  from  the  remaining  SNPs  were  used 
for  genetic  quality  control.  One  subject  from  each  pair  with  pi-hat  value  >0.2  was  randomly 
removed  (n=563),  and  outliers  from  the  principal  component  analysis  (n=94)  were  removed.  In 
total,  1408  subjects  (n=382  asthmatics,  n=1026  non-asthmatics)  and  348,697  SNPs  were 
included  in  the  replication  stage  analysis. 

In  the  WHI  SHARe,  the  subjects  were  genotypes  on  the  Affymetrix  6.0  mapping  array, 
which  contained  probes  for  approximately  1,000,000  SNPs  in  African  American  (n=8218)  and 
Hispanic  (n=3482)  subjects.  Intensity  data  generated  for  these  probes  can  then  be  utilized  to 
perform  CNV  calling,  detailed  below. 

3.4  SNP  Genetic  Analysis 

For  genetic  analysis,  SNPs  were  coded  as  an  additive  genetic  model  (0,  1,  and  2, 
indicating  the  number  of  minor  alleles),  and  BMI  was  treated  as  a  quantitative  trait.  Asthma 
status  was  dichotomously  coded.  Population  stratification  was  inspected  through  the  genomic 
inflation  factor  (7)  in  single  SNP  analyses,  with  adjustment  for  age,  sex,  asthma  status,  and 
increasing  numbers  of  principal  components.  We  found  that  the  first  principal  component  was 
sufficient  to  account  for  population  stratification.  Quantile-quantile  plots  were  based  on  p-values 
from  the  genome-wide  single  SNP  analyses.  A  genome-wide  SNP  by  asthma  interaction  analysis 
was  then  conducted  on  MESA  subjects,  with  a  linear  model  that  included  main  effect  tenns  for 
the  SNP  and  asthma  and  an  interaction  tenn  for  SNP  x  asthma.  The  model  also  included 
covariates  to  adjust  for  age,  sex,  and  the  first  principal  component.  Strict  Bonferroni-corrected 
significance  was  defined  as  p  <  6.93 xlO'8  (0.05/721,893),  while  genome-wide  suggestive 
significance  was  defined  as  p  <  5X 10'7  for  tests  of  interaction  terms.  Three  genetic  regions  that 
each  had  more  than  three  genome- wide  suggestive  SNPs  in  high  LD  (r2>0.8)  were  identified,  and 
the  SNP  in  each  region  with  the  smallest  number  of  missing  subjects  was  selected  as  a  tag  SNP 
for  replication  in  FHS. 
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In  FHS,  a  genome- wide  single  SNP  analysis  was  first  perfonned  with  a  linear  regression 
model  adjusted  for  age,  sex,  asthma  status,  and  the  first  principal  component  to  inspect  genomic 
inflation.  For  replication  of  interactions  involving  the  tag  SNPs  in  FHS,  SNP-asthma  interactions 
for  BMI  as  the  outcome  were  analyzed  using  linear  regression  models  as  already  described. 
Interactions  reaching  nominal  significance  (p<0.05)  in  FHS  were  considered  potentially 
replicated.  The  regional  plot  for  the  replicable  region,  17q2 1 .2,  was  made  using  LocusZoom 
[15].  The  rest  of  the  six  genome-wide  suggestive  SNPs  in  17q2 1 .2  were  then  tested  for 
interaction  in  FHS.  SNP  imputation  was  performed  using  IMPUTE2  [16]  with  default  settings 
for  the  two  un-genotyped  SNPs  (rs  1260 1191  and  rsl6968877)  in  FHS.  HapMap  3  haplotypes 
with  NCBI  [National  Center  for  Biotechnology  Infonnation]  build  36  (hg  1 8)  coordinates  were 
used  as  the  reference  panel.  The  imputation  had  an  overall  concordance  rate  of  96.9%,  and  the 
info  score  was  0.988  for  rs  1260 1191  and  0.979  for  rsl6968877.  All  SNPs  were  then  examined 
and  subjected  to  meta-analysis,  conducted  using  the  R  package  “meta”  [17],  with  both  fixed 
effect  and  random  effect  models.  The  most  significant  interacting  SNP  from  the  meta-analysis, 
rs2107212,  was  selected  for  the  stratification  analysis. 

3.5  Stratified  Analysis 

The  association  between  genotypes  of  rs2107212  and  BMI  was  evaluated  in  asthmatics 
and  non-asthmatics  in  MESA  and  FHS  subjects.  A  qualitative  analysis  using  obesity  (defined  as 
BMI>30  kg/m2)  vs.  non-obesity  (BMI<30  kg/m2)  as  the  outcome  was  perfonned.  A  logistic 
regression  model  with  adjustment  for  age,  sex,  and  the  first  principal  component  was  used  to 
evaluate  the  odds  ratio  and  95%  confidence  interval  (Cl)  of  being  obese  for  each  additional  risk 
allele,  stratified  by  asthma  status,  in  MESA  and  FHS  and  in  a  pooled  meta-analysis.  Since  there 
was  no  significant  heterogeneity  between  these  studies  (Cochrane’s  Q  test  p>0.05  and  I2<10%),  a 
fixed  effect  model  was  used  in  the  meta-analysis.  We  also  evaluated  the  association  when  using 
different  BMI  cut-offs  of  28  kg/m2  and  25  kg/m2  to  define  obesity.  In  FHS,  a  subset  of  subjects 
who  participated  in  exam  2  and  had  BMI  and  asthma  infonnation  was  extracted  to  examine  the 
subjects’  BMI  change  from  1979  (exam  2)  to  2005  (exam  8).  Genotype-specific  mean  BMI  at 
both  time  points  was  calculated.  Net  BMI  change  over  26  years  was  then  compared  between 
genotypes  with  or  without  the  risk  allele,  stratified  by  asthma  status.  Welch’s  t-tests  were  applied 
when  comparing  the  difference  of  two  groups. 

3.6  CNV  Detection 

Two  independent  software  packages,  PennCNV  [18]  and  Birdsuite  [19],  were  used  for 
CNV  detection  from  the  SNP  genotyping  arrays.  In  PennCNV,  CNV  calling  was  perfonned 
under  standard  protocol  and  adjusted  for  guanine-cytosine  content.  Sample  quality  controls  were 
based  on  quality  control  metrics  generated  by  PennCNV  including  Log  R  ratio  standard 
deviation,  B  allele  frequency  drift,  wave  factor,  and  total  number  of  CNVs  called  per  person. 
CNVs  that  were  less  than  10  kb  or  contained  less  than  10  SNPs  or  resided  in  the 
centromere/telomere  regions  were  excluded.  The  Birdseye  module  from  Birdsuite  was  also 
applied  to  detect  CNVs  following  standard  pipeline.  CNV  calls  with  a  limit  of  detection  score  > 

5  were  filtered  in  for  further  quality  control  procedures  in  the  same  way  as  in  PennCNV.  The 
overlapping  CNV  calls  from  both  PennCNV  and  Birdsuite  were  retained  for  further  analysis  in 
PLINK  [13]. 
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3.7  CNV  Association 


First,  we  used  two  independent  software  packages,  PennCNV  and  Birdsuite,  to  detect 
CNVs,  respectively,  in  WFII  and  then  the  overlapped  CNV  calls  were  recorded  for  further  quality 
control  in  PLINK.  In  PLINK,  CNVs  were  divided  into  duplications  and  deletions  and  filtered  out 
if  the  length  was  over  1000  kb  or  the  frequency  was  less  than  1%.  Then  CNVs  were  grouped  into 
pools  with  overlapping  segments.  There  were  68 1  groups  of  deletions  and  452  groups  of 
duplications.  Genome-wide  CNVs  by  BMI  interaction  analysis  were  performed  using  a  logistic 
regression  model  with  adjustment  for  age,  BMI  value,  and  the  first  principal  component  looking 
for  associations  with  asthma. 

4.0  RESULTS 

4.1  Summary  of  Participants 

Summary  statistics  of  participants  from  MESA  (n=2474)  and  FITS  (n=1408)  following 
quality  control  are  shown  in  Table  1.  The  distribution  of  sex,  age,  and  BMI  was  similar  for  the 
two  datasets.  Asthma  prevalence  differed  between  the  two  studies  (10%  vs.  27%  in  MESA  and 
FHS,  respectively).  BMI  was  significantly  higher  in  asthmatics  compared  to  non-asthmatics  in 
both  studies  (Welch’s  t-tests).  As  expected,  the  prevalence  of  obesity  was  significantly  higher 
among  asthmatics  compared  to  non-asthmatics  in  both  datasets  (1  degree  of  freedom  chi-square 
test). 


Table  1.  Characteristics  of  Subjects  in  MESA  and  FHS 


MESA  (n=2474) 

FHS  (n=1408) 

Variable 

Asthma 

(n=257) 

Non-Asthma 

(n=2217) 

,  Asthma 

p-ValUe  (n=382) 

Non-Asthma 

(n=1026) 

p-value 

Age  (mean±SD) 

62.91±10.12 

61.04±9.65 

68.08±8.82 

68.77±8.83 

Male  (%) 

115  (45) 

1097  (49) 

189  (49) 

476  (46) 

BMI  (mean±SD) 

29.01±6.04 

27.59±6.04 

3.25  xlO4  29.17±5.72 

27.74±4.90 

1.75xl0-s 

Classification 

Underweight 

3 

12 

4 

10 

Normal  weight 

66 

716 

82 

298 

Overweight 

89 

910 

153 

431 

Obesity  (%) 

99  (39) 

579  (26) 

3.37xl0-5  143  (37) 

287  (28) 

7.73  xlO’4 

SD  =  standard  deviation. 


Note:  Participants  who  were  included  in  the  final  interaction  analyses  are  summarized.  In  MESA,  age  is  based  on 
the  screening  exam;  in  FHS,  age  is  based  on  the  8th  exam  with  5-year  intervals.  Classification  was  based  on  BMI  as 
follows:  underweight  (BMI<18.5);  normal  weight  (18.5<BMI<25);  overweight  (25<BMI<30);  obesity  (BMI>30). 
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4.2  Genome-Wide  SNP  by  Asthma  Interaction  Analysis 

Genome-wide  SNP  by  asthma  interaction  analysis  was  first  performed  on  MESA 
subjects.  The  linear  regression  model  included  main  effects  for  SNP  and  asthma  and  an 
interaction  term  for  SNP  x  asthma.  The  model  was  also  adjusted  for  age,  sex,  and  the  first 
principal  component.  Overall  interaction  results  are  shown  in  the  Manhattan  plot  (Figure  1). 
Thirty-one  interacting  SNPs  exceeded  the  genome- wide  suggestive  threshold  of  p<5><10"7,  and 
three  genetic  regions  that  each  contained  more  than  three  genome -wide  suggestive  interacting 
SNPs  in  high  LD  with  each  other  were  selected:  7p21.3,  10q25.3,  and  1 7q2 1 .2.  The  SNPs  with 
the  least  missing  subjects  in  each  region  were  chosen  as  tag  SNPs  for  further  validation  in  FHS. 
Of  the  three  tag  SNPs,  one  SNP  (rs2107212  in  1 7q2 1 .2)  had  a  significant  SNP  by  interaction 
p-value  in  replication  subjects  (p=9.38x  10'3),  while  the  other  two  SNPs  did  not  reach  the 
replication  significance  threshold  of  0.05  (Table  2).  To  further  verify  the  significant  region,  we 
then  tested  all  seven  MESA  genome-wide  suggestive  SNPs  in  1 7q2 1 .2  in  FITS.  We  observed  that 
all  seven  SNPs  had  an  SNP  by  asthma  interaction  replication  p-value  <  0.05,  in  the  same 
direction  as  observed  in  MESA.  The  p-values  for  the  seven  SNPs  ranged  from  1 .76x  10"7  to 
3.72x  10"7  in  MESA  and  9.38x  10"3  to  0.047  in  FHS.  The  regional  plot  for  SNPs  located  within 
200  kb  of  the  tag  SNP,  rs2107212,  in  MESA  is  presented  (Figure  2),  where  a  clear  cluster  of 
SNPs  in  strong  LD  is  shown.  SNP  rs2107212  was  also  the  most  significant  asthma-interacting 
SNP  in  this  region  following  a  meta-analysis  of  the  MESA  and  FHS  data  (p=2.5 1  x  10"3  under  a 
random  effect  model  and  p=5.60x  10'8  under  a  fixed  effect  model). 


Chromosome 


Figure  1.  Manhattan  plot  of  interaction  p-values  derived  from  the  MESA  genome-wide  interaction  analysis. 

Linear  regression  models,  including  main  effects  for  SNP  and  asthma  and  an  interaction  term  for  SNP  x  asthma, 
were  used.  Models  were  also  adjusted  for  age,  sex,  and  the  first  principal  component.  The  red  (upper)  line 
represents  the  genome-wide  significant p-  value  (0.05/721, 893=6.93 x10's).  The  blue  (lower)  line  represents  the 
genome-wide  suggestive  p-value  of  5  xlO'7. 
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Table  2.  Top  Interacting  Regions  Associated  with  BMI  and  Replication  in  FHS 


SNP 

Region 

Position 

Gene/Location 

Dataset 

p  (95%  Cl) 

p-value 

rs  10250689 

7p21.3 

10120482 

HSPA  8P8/3  3  Okb  up 

MESA 

4.66  (2.98,  6.34) 

5.91E-08 

FHS 

0.56  (-0.98,  2.10) 

4.75E-01 

rs7904383 

10q25.3 

115716396 

ATRNLl/Intron 

MESA 

-2.66  (-3.63,  - 
1.70) 

7.42E-08 

FHS 

-0.06  (-0.94,  0.81) 

8.87E-01 

rs2107212 

1 7 q2 1 . 2 

40949109 

KRT39/9kb  down 

MESA 

2.96(1.82,4.09) 

3.72E-07 

FHS 

1.49  (0.37,2.61) 

9.38E-03 

P  =  interaction  coefficient. 

Note:  From  the  MESA  gene  by  environment  interaction  analysis,  regions  with  more  than  three  genome-wide 
suggestive  SNPs  (interacting  p-value  <5X10"7)  in  strong  LD  (r2>0.8)  were  selected  as  top  interacting  regions.  SNPs 
with  fewer  missing  subjects  in  each  region  were  chosen  as  tag  SNPs  to  replicate  in  FFIS.  SNP  position  was  based  on 
the  GRCh38  assembly. 
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Figure  2.  Regional  plot  showing  replicable  SNPs  in  17q21.2.  A  regional  plot  was  generated  in  LocusZoom  v.  1.1 
for  the  MESA  genome-wide  interaction  analysis  results,  p-values  on  the  -loglO  scale  are  displayed  on  they  axis, 
and  chromosome  positions  are  displayed  on  the  x  axis.  This  plot  shows  the  tag  SNP  (rs2107212)  in  purple  with  a 
200-kb  flanking  region  on  each  side.  The  pairwise  LD  pattern  with  rs2 107212  is  also  shown.  The  bottom  panel 
shows  the  name,  position,  and  transcription  direction  of  each  gene  in  this  region. 
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4.3  Association  Between  BMI  and  rs2107212  Genotypes 


To  further  explore  the  nature  of  the  interaction,  we  conducted  a  series  of  stratified 
analyses  for  rs2107212  among  the  asthmatic  and  non-asthmatic  populations  in  MESA  and  FHS. 
First,  we  evaluated  the  association  between  rs2107212  and  BMI  stratified  by  asthma  status  in 
MESA  (Figure  3A).  Among  asthmatics,  increasing  copies  of  the  A  allele  at  rs2107212  were 
associated  with  successively  higher  BMI  values.  Interestingly,  among  non-asthmatics,  an 
opposite  trend  was  observed,  although  the  differences  were  not  significant.  Subjects  carrying  the 
AA  or  AG  genotype  had  significantly  higher  BMI  among  asthmatics  compared  to  non¬ 
asthmatics.  In  FHS,  we  observed  a  similar  trend  for  the  effect  of  the  A  allele  (Figure  3B). 
Subjects  carrying  the  AG  genotype  had  a  significantly  greater  BMI  difference  between 
asthmatics  and  non-asthmatics  compared  to  those  carrying  the  GG  genotype.  Subjects  carrying 
the  AA  genotype  had  higher  BMI  among  asthmatics  compared  to  non-asthmatics,  although  the 
p-value  did  not  reach  nominal  significance  (p=0.06). 


A 


8  n 


non- asthmatic?  asthmatics 


□  AA 

□  AG 

□  GG 


group 
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8  -i 

S  - 
8  - 
8  - 


□  AA 

□  AO 

□  00 


group 

Figure  3.  Association  of  rs2107212  genotypes  and  BMI,  stratified  by  asthma  status  in  MESA  (A)  and  FHS  (B), 
respectively.  BMI  (kg/m2)  was  compared  across  subjects  carrying  the  AA,  AG,  or  GG  genotypes  for  rs2107212  in 
asthmatics  and  non-astlvnatics.  Numbers  of  subjects  are  shown  in  each  bar.  Data  are  represented  as  the  mean  ± 
standard  error  measurement.  Significance  was  tested  using  Welch ’s  t-test  (*p<0.05,  **p<0.01,  and  ***p<0.001). 
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4.4  Evaluation  of  the  Odds  Ratio  of  Being  Obese  for  the  A  Allele  at  rs2107212 

We  next  examined  and  compared  the  odds  ratio  (OR)  of  being  obese  (BMI>30  kg/nr) 
among  the  asthmatic  and  non-asthmatic  populations  for  rs2107212  in  MESA  and  FHS  separately 
as  well  as  together  in  a  meta-analysis.  The  forest  plot  showed  a  consistent  effect  of  the  A  allele 
in  MESA  and  FHS  (Figure  4).  The  meta-analysis  showed  that  the  overall  odds  of  being  obese 
increased  by  1.89  fold  for  each  additional  A  allele  in  the  asthmatic  population.  In  the  non¬ 
asthmatic  population,  the  odds  of  being  obese  decreased  by  0.89  fold  per  additional  A  allele; 
however,  this  was  not  significantly  different  from  1.0.  To  further  test  if  the  asthma-specific 
association  results  for  the  A  allele  at  rs2107212  are  still  present  with  varying  cut-offs  for 
classifying  obesity,  we  evaluated  the  association  using  different  BMI  cut-offs  of  28  kg/m2  and 
25  kg/m2.  These  results  showed  that  different  BMI  cut-offs  resulted  in  similar  trends  of 
association  for  rs2107212.  This  further  demonstrated  that  the  A  allele  at  rs2107212  was 
associated  with  BMI  increase  and  obesity  incidence  only  among  asthmatics. 


Asthmatic  group 
Study  OR(95%CI) 
MESA  2.30(1.43,3.69) 
FHS  1.65(1.11.2.46) 
Meta  1.89(1.39,2.57) 
Heterogeneity: 
p=0.30  l-squared=7.3% 


Non-asthmatic  group 
Study  OR(95°/.Cl) 
MESA  0.87(0.73.1.04) 
FHS  0.94(0.73.1.21) 
Meta  0.89(0.77,1.03) 
Heterogeneity: 
p=0.62  l-squared=0% 
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Figure  4.  Forest  plot  showing  the  association  of  rs2107212  with  obesity  in  MESA,  FHS,  and  meta-analysis. 

Asthma  status-specific  ORs  of  being  obese  for  rs2107212  were  calculated  using  logistic  regression  models  adjusted 
for  age,  gender,  and  the  first  principal  component  in  MESA  and  FHS,  respectively.  The  minor  allele  A  is  the  “risk” 
allele,  and  an  additive  genetic  model  was  used.  Meta-analysis  was  performed  under  a  fixed  effect  model  since  the 
heterogeneity  tests  were  not  significant.  Boxes  indicate  group-specific  OR  point  estimates,  and  lines  indicate  the 
respective  95%  CL  Diamonds  indicate  meta-analysis  OR  and  95%  Cl. 

4.5  BMI  Change  Over  a  26-Year  Period  After  Asthma  Diagnosis  in  FHS 

In  FHS,  990  subjects  had  BMI  and  asthma  data  for  both  exams  2  (1979)  and  8  (2005),  of 
which  983  subjects  had  genotype  information  for  rs2107212.  Mean  BMI  at  these  two  times, 
stratified  by  asthma  status  and  rs2107212  genotype,  is  shown  in  Figure  5A.  BMI  among 
asthmatics  was  significantly  higher  than  in  non-asthmatics  in  1979  (p=0.0102),  and  the 
difference  was  more  significant  in  2005  (p=6.84x  1 0"4).  During  this  26-year  period,  in  asthmatics, 
subjects  carrying  the  AA  or  AG  genotype  increased  their  BMI  more  than  subjects  carrying  the 
GG  genotype;  however,  this  was  not  true  among  non-asthmatics  (Figure  5B).  During  this  time 
period,  asthmatics  carrying  the  AA  or  AG  genotypes  had  a  mean  BMI  change  of  3.79,  while 
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subjects  carrying  the  GG  genotype  had  a  mean  BMI  change  of  3.24.  Among  non-asthmatics, 
subjects  carrying  the  AA  or  AG  genotypes  had  a  mean  BMI  change  of  2.3 1,  while  subjects 
carrying  the  GG  genotype  had  a  mean  BMI  change  of  2.71.  This  result  further  indicates  that  the 
asthmatic  population  in  general  gained  more  weight  than  the  non-asthmatic  population  and  that 
the  A  allele  at  rs2107212  increased  the  BMI  change  by  60%  among  asthmatics  compared  to  non¬ 
asthmatics  (p=0.0308). 


A 


B 


year 


■  AA  AG 
□  GG 


group 


Figure  5.  Association  of  rs2107212  genotypes  and  BMI  change  by  asthma  status  in  FHS.  (A)  For  a  subset  of 
FHS  subjects  who  had  BMI  and  asthma  information  at  both  exam  2  (1979)  and  exam  8  (2005)  (n=983),  the  means 
of  BMI  at  both  time  points  are  shown  according  to  the  rs2107212  genotypes  AA/AG  and  GG,  stratified  by  asthma 
status.  (B)  Net  BMI  change  (BMI  value  in  2005  —  BMI  value  in  1979)  with  rs2107212  genotypes  AA/AG  and  GG 
stratified  by  asthma  status  is  shown.  The  number  of  subjects  is  shown  in  each  bar.  Data  are  represented  as  the  mean 
±  standard  error  measurement.  *p<0.05  under  Welch ’s  t-test. 


4.6  CNV  Results 

We  have  identified  top  CNVs  for  African  Americans  in  the  WHI  dataset.  The  most 
significant  deletion  interacting  with  BMI  for  asthma  is  on  chromosome  7,  from  1 17927028  base 
pairs  (bp)  to  117944067  bp  (p-value  =  0.0032;  Table  3).  The  most  significant  duplication  region 
that  is  interacting  with  BMI  in  the  outcome  of  asthma  is  in  chromosome  10,  from  46496693  bp 
to  46499590  bp  (p-value  =  0.0019;  Table  4). 
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Table  3.  Top  5  Deletions  in  Interaction  Analysis 


Chromosome 

Starting  bp 

Ending  bp 

OR 

95%  Cl 

p-value 

7 

117927028 

117944067 

1.07 

1.02,  1.12 

0.0032 

1 

110044476 

110044476 

1.18 

1.04,  1.33 

0.0080 

9 

43666998 

43666998 

1.08 

1.02,  1.15 

0.0119 

2 

35437965 

35459760 

1.07 

1.01,  1.13 

0.0130 

1 

110038455 

110038470 

1.15 

1.03,  1.28 

0.0146 

Table  4.  Top  5  Duplications  in  Interaction  Analysis 


Chromosome 

Starting  bp 

Ending  bp 

OR 

95%  Cl 

p-value 

10 

46496693 

46499590 

1.11 

1.04,  1.19 

0.0019 

10 

46501701 

46501866 

1.11 

1.04,  1.19 

0.0024 

10 

46504791 

46504871 

1.11 

1.04,  1.18 

0.0030 

10 

46506308 

46506308 

1.11 

1.04,  1.18 

0.0030 

10 

46510869 

46510869 

1.10 

1.03,  1.17 

0.0039 

5.0  DISCUSSION 

Many  correlated  disorders  have  shared  genetic  backgrounds,  and  the  genetic  effects  on 
one  of  the  phenotypes  might  be  modified  by  the  other  phenotype.  This  provides  a  novel  model 
supporting  the  integration  of  gene  by  environment  interaction  terms  (i.e.,  where  the 
“environment”  is  another  disease  or  phenotype)  in  a  GWAS  framework  to  explore  comorbid 
susceptibility  genes.  These  genetic  factors  could  be  overlooked  when  only  SNP  main  effects  are 
examined. 

The  comorbidity  of  asthma  and  obesity  is  a  growing  medical  problem.  While  previous 
studies  have  focused  on  the  unidirectional  relationship  of  obesity  preceding  asthma  onset,  post¬ 
asthma  BMI  increase  is  a  neglected  but  potentially  important  phenomenon,  particularly  since 
asthma  is  an  early  onset  disease.  Here,  we  used  two  independent  datasets  of  older  adults  and 
hypothesized  that  asthma  impacted  BMI  through  genetic  factors.  Therefore,  we  conducted  an 
SNP  by  asthma  interaction  analysis  to  identify  SNPs  that  are  associated  with  BMI  and  modified 
by  asthma  status.  SNPs  in  1 7q2 1 .2  were  verified  to  be  associated  with  higher  BMI  and  higher 
risk  of  obesity  only  when  asthma  is  present. 

We  further  took  advantage  of  the  longitudinal  FHS  dataset  to  obtain  BMI  measured 
closest  to  asthma  onset  and  evaluated  BMI  change  after  asthma  diagnosis  from  1979  to  2005. 
Figure  5  A  reflects  the  trend  of  BMI  change  over  this  26-year  period.  In  1979,  the  asthmatic 
population  already  had  significantly  higher  BMI  than  non- asthmatics.  The  BMI  differential 
further  increased  in  the  intervening  26  years  and  resulted  in  a  much  greater  BMI  differential 
among  asthmatics  than  non-asthmatics.  When  genotype  was  taken  into  consideration,  asthmatic 
subjects  carrying  the  minor  allele  of  rs2107212  gained  more  weight  than  subjects  not  carrying 
this  allele,  while  this  was  not  the  case  among  non-asthmatics.  Since  BMI  in  1979  was  likely 
measured  after  but  closer  to  asthma  onset  than  BMI  in  2005,  we  can  infer  that  the  BMI 
difference  between  asthmatics  and  non- asthmatics  before  asthma  onset  was  less  than  in  1979  and 
that  actual  BMI  increase  after  asthma  may  be  greater.  In  Figure  5B,  we  further  examined  the  net 
BMI  change  over  a  26-year  period,  where  rs2107212  functioned  as  a  risk  locus  for  BMI  increase 
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only  when  asthma  was  present.  Here,  we  identified  a  genetic  locus  that  is  associated  with  post¬ 
asthma  BMI  increase  and  demonstrated  that  asthma  could  influence  weight  gain  over  a  long 
period. 

Chromosome  1 7q2 1  has  been  repeatedly  reported  to  be  associated  with  asthma  [20-22], 
Although  the  region  covering  the  ORMDL3,  GSDMA,  and  GSDMB  genes  is  the  most  replicated 
one,  the  keratin  (KRT)  cluster  has  also  been  shown  to  be  associated  with  asthma  in  a  previous 
GWAS  study  [23],  KRT  proteins  are  key  structural  components  in  epithelial  cells  and  are 
essential  for  tissue  function.  A  study  identified  KRT  1 8  as  a  bronchial  epithelial  autoantigen  that 
is  associated  with  adulthood  nonallergic  asthma  [24].  KRT  18  and  some  other  keratin  genes  have 
also  been  identified  to  be  associated  with  obesity-caused  fatty  liver  [25,26].  Hair  follicles  have 
been  demonstrated  to  be  associated  with  Toluene  diisocyanate-induced  asthma  [27].  Non¬ 
functional  mutations  in  filaggrin,  a  protein  binding  to  and  condensing  the  keratin  cytoskeleton, 
are  risk  factors  for  asthma,  atopic  eczema,  and  allergies  [28].  The  identified  SNP  rs2107212  is 
located  in  the  intergenic  region  flanked  by  keratin  genes  KRT39  and  KRT23  (9  kb  upstream  of 
KRT39  and  1 1  kb  downstream  of  KRT23).  KRT39  is  a  type  I  keratin  gene  with  expression  late 
in  the  differentiation  of  hair  [29].  KRT23  is  a  type  I  epithelial  keratin  gene  that  has  been 
identified  to  interact  with  14-3-3  proteins  to  modulate  key  cellular  processes  in  a  SMAD4- 
dependent  manner  [30].  In  addition,  KRT23  has  been  reported  to  be  associated  with  several 
cancers,  such  as  colon  cancer,  pancreatic  cancer,  and  hepatocellular  carcinoma  [31-33].  This 
study  links  the  keratin  family  genes  with  asthma  and  obesity  phenotypes.  However,  the 
underlying  biological  mechanisms  of  rs2107212  in  asthma-dependent  BMI  increase  need  to  be 
further  explored. 

This  study  is  potentially  confounded  by  reverse  causality  from  obesity  to  asthma. 
Longitudinal  studies  have  observed  an  increased  OR  or  relative  risk  of  incident  asthma  among 
the  obese  population  [34,35].  We  cannot  exclude  the  possibility  of  a  causal  effect  of  obesity  on 
asthma,  since  BMI  before  asthma  development  was  not  available  in  MESA  and  the  Framingham 
offspring  study.  However,  to  attempt  to  further  assess  SNP  effects  on  BMI  increase  after  asthma 
diagnosis,  we  used  a  subset  of  subjects  from  the  Framingham  offspring  cohort  who  had  both 
asthma  status  and  BMI  values  in  1979  and  2005.  Over  a  period  of  26  years,  the  net  BMI  change 
was  significantly  higher  in  asthmatics  than  non-asthmatics.  This  result  demonstrates  that  SNP 
rs2107212  is  involved  in  the  effect  of  asthma  on  BMI.  However,  from  the  data  we  have,  we 
cannot  detennine  if  rs2107212  is  also  involved  in  the  reverse  relationship,  BMI  on  asthma. 

Another  limitation  in  this  study  is  that  MESA  and  FHS  did  not  use  the  same  asthma 
definitions.  MESA  asked  about  doctor-diagnosed  asthma,  while  FHS  asked  about  wheezing  and 
asthma  in  the  same  question  and  a  doctor  diagnosis  was  not  required.  This  could  explain  the 
asthma  prevalence  difference  (10%  in  MESA  and  27%  in  FHS).  A  relatively  looser  asthma 
definition  in  FHS  may  have  potentially  diluted  any  signals  and  may  have  underestimated  some  of 
the  interactions,  which  we  believed  was  the  primary  reason  for  lack  of  replication  of  other  top 
interacting  regions.  However,  the  consistency  in  direction  of  effect  between  MESA  and  FHS  for 
the  tag  SNPs  strongly  argues  in  favor  of  these  being  true  positive  interactions. 

After  adjusting  for  the  number  of  SNPs  tested  in  the  genome-wide  screen,  this  study  had 
over  80%  power  to  detect  a  P  of  5  when  the  SNP  MAF  is  >0. 1  and  a  P  of  4  when  the  SNP  MAF 
is  >0.2.  While  this  is  a  modest  level  of  power,  this  dataset  allowed  us  to  identify  a  novel  and 
replicable  gene  x  environment  interaction. 
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The  results  of  the  CNV  by  BMI  analysis  in  African  Americans  in  the  WHI  study  were  not 
genome-wide  significant  when  looking  for  associations  with  asthma.  However,  the  top  deletion 
and  duplication  were  nominally  p<0.05  and  may  be  worth  further  exploration  in  the  future  in 
other  datasets  containing  African  Americans  and/or  other  ethnic  groups. 

6.0  CONCLUSIONS 

In  conclusion,  this  is  the  first  attempt  at  a  genome-wide  interaction  analysis  to  detect 
genetic  risk  factors  associated  with  BMI  modified  by  asthma.  SNPs  in  17q21.2  were  identified  as 
risk  loci  for  BMI  increase  only  among  asthmatics.  This  finding  will  help  elucidate  pathways 
involved  in  the  comorbidity  of  asthma  and  obesity.  Further,  we  identified  some  candidate  CNVs 
that  interact  with  BMI  that  are  associated  with  asthma.  Additional  studies  are  needed  to  further 
verify  and  emphasize  the  relationship  between  asthma,  obesity,  and  genetic  markers  and  examine 
the  underlying  mechanisms. 
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LIST  OF  ABBREVIATIONS  AND  ACRONYMS 


p 

interaction  coefficient 

BMI 

body  mass  index 

bp 

base  pairs 

CDI 

clinical  diagnostic  impression 

Cl 

confidence  interval 

CNV 

copy  number  variant 

FHS 

Framingham  Heart  Study 

GWAS 

genome-wide  analysis  study 

KRT 

keratin 

LD 

linkage  disequilibrium 

MAF 

minor  allele  frequency 

MESA 

Multi-Ethnic  Study  of  Atherosclerosis 

OR 

odds  ratio 

SHARe 

SNP  Health  Association  Resource 

SNP 

single  nucleotide  polymorphism 

WHI 

Women’s  Health  Initiative 
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